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A truncated partial wave analysis for the photoproduction of two pseudoscalar mesons on a nucleon 
is discussed with respect to the determination of a complete set of observables. For the selection of 
such a set we have applied a criterion previously developed for photo- and electrodisintegration of 
a deuteron, which allows one to find a 'minimal' set of observables for determining the partial wave 
amplitudes up to possible discrete ambiguities. The question of resolving the remaining ambiguities 
• by invoking additional observables is discussed for the simplest case, when the partial wave expansion 

' is truncated at the lowest total angular momentum of the final state Jmax = 1/2. The resulting 

'fully' complete set, allowing an unambiguous determination of the partial wave amplitudes, is 
£\J ■ presented. 

$—1 '. 

PACS numbers: 13.60.Le, 13.75.-n, 21.45.+V, 25.20.Lj 

V/-) \ I. INTRODUCTION 

(N ; 

Photoproduction of two pseudoscalars on nucleons have been studied rather intensively during the last two decades. 
At present a large amount of experimental data, in particular for 7r°7r° and Tr°r] channels, have been collected and 
some new experiments on polarization observables for these reactions are planned. The interpretation of the data 
within different models has allowed one to qualitatively understand the major mechanisms of these processes. At the 
same time, although the general agreement of the various calculations with the measured cross sections is reasonable, 
significant qualitative differences between these models exist. Some of them were already discussed, for example in 
1 Ref. Q. 

One reason for these differences between theoretical results is that the standard approach, based on the isobar 
model, appears to have reached certain limitations. Probably its weakest point is that the corresponding formalism 
I , depends on a specific assumption about the dynamics of the production process. Within the isobar model approach 
one assumes that the two-body discontinuity in the reaction matrix coming from the interaction in the two-body 
subsystems in the final state may be approximated by resonance terms (usually taken in the Breit-Wigner form). In 

■ other words, it is assumed that the final three-particle state is produced via intermediate formation of quasi two-body 
£f~) [ states containing meson-nucleon and meson-meson isobars. 

In order to achieve a significant improvement of present theoretical models one has to eliminate as much as possible 
the mentioned model dependencies from the formal description of these reactions. In this respect, an ideal tool for 
the investigation of the reaction dynamics is an analysis of a complete experiment within a given model, based on the 

■ fundamental principles of rotation and parity invariance. It is clear that for the photoproduction of two mesons this 
task is considerably more complicated in comparison to the photoproduction of a single meson. Firstly, in the case of 
three particles in the final state, the amplitudes depend on five kinematical variables, so that for their determination 
accurate measurements of five-dimensional distributions are needed. Secondly, contrary to the single meson case, a 
complete set contains a considerably larger number of observables. Indeed, if two mesons are produced, the property 
of parity conservation does not allow one to reduce the number of independent amplitudes. Therefore, in order to 
determine all eight complex amplitudes up to an overall phase one needs at least 15 observables. As was shown by 
Roberts and Oed Q, a complete set includes not only single and double but also triple polarization observables. The 
latter is especially disappointing since it requires the implementation of complicated measurements at a high level of 
accuracy. 

At the same time, as was noted in Refs. Q and [3], arbitrariness in the overall phase at each point of the phase space 
does not allow one to find the multipole amplitudes, which are obviously needed for a nucleon resonance analysis. 
Therefore, if one searches for resonances or, more generally, for states with definite spin and parity J p , it is more 
reasonable (and probably less complicated) to adopt a truncated partial wave expansion up to a maximal total angular 
momentum J max and to study instead of the spin or helicity amplitudes the partial wave amplitudes. In this case, as 
a rule, a lower number of polarization observables is needed. As is discussed in Refs. for 'fN — > n°N, in order 

to determine (up to an overall phase) the multipoles Eq + , E\ + , and M%± which are important in the first resonance 
region, already a set of single polarization observables with an additional measurement of only one of the double 
polarization observables (for example, F- or G-asymmetry) is sufficient. 
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For a partial wave analysis we firstly need a convenient and model independent form of the partial wave decompo- 
sition of the reaction amplitude. As already noted above, the isobar model does not meet this requirement, since it 
depends on specific assumptions about the reaction dynamics. Therefore, we adopt in the present paper the formalism 
developed in Ref. 0- Here a special coordinate frame is used in which the z-axis is chosen along the normal to the 
plane spanned by the momenta of the final particles in the overall center-of-mass frame. In this case one can choose 
as independent kinematical variables the energies of two of the three final particles and three angles, determining the 
orientation of the final state plane with respect to the photon beam (one of the angles corresponds to the rotation 
around the normal). Then the amplitude for a given total angular momentum J of the final state and its projection 
M on the z-axis is obtained by an expansion of the helicity amplitudes over a set of Wigner functions. In Ref. [lj this 
approach was applied successfully to the analysis of the unpolarized differential cross section for the photoproduction 
of 7T 7r° pairs. 

The second important question is what is an optimal choice of the obscrvablcs needed for the determination of 
the partial wave amplitudes. While for the determination of n complex quantities (up to an overall phase factor) 
only In — 1 real parameters are needed, the total number of linearly independent observables (or, more generally, the 
number of bilinear hermitian forms of the amplitudes) is n 2 . This means that there exist also nonlinear dependencies 
between the observables, so that not any subset of 2n — 1 observables may form a complete set from which the In — 1 
real parameters can be deduced. 

Thus the question is how to select out of this set of n 2 observables a subset of 2n — 1 independent observables. In 
the present paper we describe a method which may be used to single out a minimal set of independent observables. 
The criterion, which underlies the method, was originally developed for deuteron photo- and electrodisintegration in 
Refs. [8, 9]. An additional very important question concerns possible discrete ambiguities, which naturally can appear 
since the extraction of the amplitudes implies a solution of a set of quadratic equations. To resolve these ambiguities 
we adopt for our truncated partial wave analysis an additional criterion, which was used in Ref. [lOj for the spin 
amplitude analysis of single meson photoproduction. 

First, we present in the next section a brief review of the formalism developed in Ref. @. Then we give in Sec. IIIII 
an example of the method of how to find a complete set of observables for the simplest case, when only partial waves 
with J = 1/2 and both parities are included. In Sect. IIVI we summarize our results and give an outlook on future 
developments. Some details are collected in two appendices. 

II. FORMAL DEVELOPMENTS 

For the theoretical description we choose the overall cm. system. The four-momenta of incoming photon, outgoing 
mesons, and final nucleon are denoted by (w 7 , k ), (wx, qi ), (u>2, 92 ), and (E 7 p), respectively. We consider within this 
system two right-handed orthogonal coordinate systems: (i) one associated with the incoming photon called with 
z-axis along the photon momentum and x-axis arbitrary, and (ii) the so-called "rigid body" system Kf Sl associated 
with the final state plane spanned by the final three particles, in which the z-axis is taken to be the normal to this 
plane and parallel topx q\. Thus the x- and y-axes are in the final scattering plane (see Fig. [l}. The transformation 
from Kf s to K 1 is given by a rotation through Euler angles (<p,9,0). Thus relative to Kf s the photon momentum k 
has the spherical angles fi 7 = (0 7 , 7 ). 

For the T-matrix we had derived in Ref. Q the following expression expanding the final state into partial waves 

r^o^^rg = e-^Y, t i^ J (^^)D( IjX _^R 1P ), (i) 

JM., 

where J and Mj denote respectively the total angular momentum of the partial wave and its projection on the 
normal to the final state plane. The rotation matrix D J M , M is taken in the convention of Rose with argument 
R-yp = (0 7P , Oj, —4>j) and <p lp = 7 — <j) p . The helicities of photon and initial and final nucleons are denoted by A, u, 
and v, respectively. As independent variables we had chosen besides the photon angles f2 7 = (0 7 , 7 ) and the proton 
angle <j) p in the final state plane, the energies of the two final mesons u)\ and UJ2- The latter two determine the relative 
angle cj> qp = <p q — cj> p between the momentum of the final proton and the relative momentum q = (q\ — 9*2 )/2 of the 
two mesons according to 

cos 4>„ = J- (lo 2 -co 2 - M 2 + M 2 ) , (2) 

with the final nucleon momentum 



P=W\ = ^{W -^-Lazf-M 2 



(3) 
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FIG. 1: Definition of the coordinate system in the cm. system. 



where W denotes the invariant total energy, and the relative momentum q of the two mesons is determined by 

q 2 = \{ul+ul-M*-Mi)-^. (4) 

In the foregoing equations the masses of the nucleon and the two mesons are denoted by Mjv, Mi and M2, respectively. 

The expression in Eq. ([T]) is obtained by making use of rotation and inversion invariance (that is angular mo- 
mentum and parity conservation) and is therefore completely general. The final partial wave is taken in the form 
\qp; ((l p ^)jpl q )JMj)(~' > where l p and l q denote the angular momenta of nucleon and meson pair, respectively. 

The contribution of the final partial wave to the reaction amplitude is given by 

Oi, W2 ) = o<M 

£ ( I -v ) ( p. - A A X ) d -«>/ 2 ) e^'- m ^ O x M L f(l p j p m p ) , (5) 



with 



X M L f(l P j P m p ) = ( ^ J E i L (-l) l * +j * +l «Tp7 P T q L$ mg (n/2) 

l q m q 



(i i .QiPT^iMn^), (6) 



where I = + 1 and O xl denotes the electromagnetic multipole operator with electric and magnetic contributions 

Om = E L M + \Mm . (7) 

The amplitudes t^f, -7 obey the following symmetry property which follows from parity conservation (see Ref. Q, 
one should note a misprint in the phase) 

tL?l x ^(K) = (-i)"- Mj *^ Mj (-^) ■ (8) 

The interchange 4> qP —<l>qp corresponds to the interchange (iVi,Mi) O (w 2 , M 2 ), see Eq. @ through Qj. Therefore, 
parity conservation does not reduce the number of independent amplitudes in this case in contrast to single meson 
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production as already noted in Ref. Q. However, in a more general sense this relation allows one to reduce the number 

JMj 
vA/i 



of independent amplitudes. Namely, provided that the amplitude t v ^ is known in the whole region of the Dalitz plot 



(tjj\,LJ2) the amplitude t J _ v _£_ can be obtained using the symmetry relation of Eq. ([5]). Thus in this more general 
sense one has 4(2J + 1) independent amplitudes for a given J, except for J = 1/2 for which this number is reduced 
to 2(2,7 + 1) = 4 because of an additional requirement A — /i < </, coming from angular momentum conservation. 

As shown in Appendix A, one can separate the contributions of those final states with positive parity from those 
of negative one according to 

C J = H*-* ( * +M) (C J + + 2 < Mj ") ■ 0) 

where 

An interesting consequence of Eq. ([9]) is that the bilinear expression t^x/a^txa with v ' — v is invariant under the 
transformation (parity exchange) 

which, as follows from Eq. (fTQ)) . is equivalent to the transformation of the matrices t J J£j 

tlZ' -> iff. (12) 

Indeed, one finds successively (note 2v = ±1) 

..I'M',* JMj 



= &u 


2uf J ' M / 








= {to4f- 




T(2"C" 


f- , j. JMj 
' c A/j 


") 




^ l AV 




- , *«J 
~r ^A/i 


+ \ _ t J'M'j*.JMj 



(13) 

Therefore, in order to distinguish between the contributions of states with different parities one has to measure recoil 
polarization along the x or y axes which are governed by such bilinear expressions with v 1 ^ v. 

As observables we will consider the differential cross section and the recoil polarization for unpolarized and circularly 
polarized photons. As shown in Ref. [l2l | the differential cross section with circular beam asymmetry is given by 



d 5 a(P?) _ da, 



l + i?E c ), (14) 



d(f)pdujiduj2dQ, 1 d(j) p duiiduj2d£l 1 
where denotes the degree of circular polarization. The unpolarized differential cross section is 

d(T ° = T o ^ ( 15 ) 

d(j) p dujiduj2di}-y 00 ' 

and the beam asymmetry for circular photon polarization 

£ C T ° = T C , (16) 

where 

TV C = c(W)Rev° £ 00 (17) 

with c(W) = M 2 N /{A{2ir) A (W 2 - M 2 N )) as a kinematical factor. Here the quantities v< f' m' -im with (x €E {0,c} are 
defined by 



v° - 
V I'M';IM — 

■ '.wo A 



V I'M':IM -] I r 

1 + OA/0 ' , 



rrT-E^WM. ( 18 ) 

1 + OMO ^ 

k~ £ Au WAf. ( 19 ) 
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where the Uj,^,. IM contain bilinear combinations of the partial wave amplitudes t J v ^ , given in Eq. ([5]), according to 



7/ VA 

U I'M';IM 



(-) M I'T^2(2j + l)Di n , m (Il 1P ) 



X 



x £ ( _ ir M. ; f /' / ./ ^ / •/' / / 

J'M',JMj 



M'j -Mj m'l I A' - m' M - A m 



x^,^,^)*^^,^), (20) 
with i? 7)9 = (</> 7P , # 7 , —<fij)- These quantities have the following symmetry property: 

4 X M<;IM = (-1) M ' +M (4-M<;I-m)* • (21) 

One should note that in Eqs. (fT5|) and (TK)|) only u^,^ t ,. IM with A' = A appear. For the recoil polarization component 
P Xi of the outgoing nucleon one has 



p d 5 a(P?) _ da 



u,4>pdcJidcJ2dO,r 1 d(j) p duiiduj2dQy 
with recoil polarizations for unpolarized beam and target 



(p° + (22) 



r xi J 00 — - r 00 ) y^ ) 
as well as beam asymmetries for circularly polarized photons 

P X C < T ° = P ^ C , (24) 

where for a € {0, c} 

(P/Q)So" - -VicWRe/Im^+o, (25) 

(W)fc a = TV2c(^)Im/Re<i- 00 , (26) 

(P/Q)ocf = c(W)Re/lmw^+ , (27) 

with 



"I'M'-JM 



\ (e iM '*» vf, M , IM ± R M ' e- iM '^ «?_ MV m) • (28) 



This completes the formal part. 



III. A TRUNCATED PARTIAL WAVE ANALYSIS 



In this section we consider a method which allows one to find for a reaction with n independent complex amplitudes 
a complete subset of 2n — 1 independent observables. It was developed in Refs. p| Q and applied to the analysis of 
deuteron electro- and photodisintegration, and we refer the reader to this paper for more details. At first, in order to 
explain the key points of the method, we consider a very simple mathematical example. Given a 2-dimensional real 
vector x — {xi, X2}, whose components are called "amplitudes" , and two real quadratic forms f^, called "observables" , 

1 2 

A(f) = -^i4i j) fc = l,2, (29) 

ij=l 

where the two matrices A k are symmetric, then the question is under which conditions for the matrices one can 
determine the amplitudes x 2 } from given values {fx, f 2 } of the observables. In other words, what is the criterion, 
that the set of quadratic equations (|29[) can be inverted (apart from possible quadratic ambiguities). 
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For the moment being, let us assume that x° = {x\,x 2 } is the required solution of Eq. (|29|) for given values 
{/i>/2j- A necessary condition for the inversion is that in the neighborhood of x° the Jacobian of the transition 
{xx, x 2 } — > {fx, h} is nonvanishing, i.e. using Eq. (|29|) 

det ( 21* ) = det 

\ dx i/ \ A\ x xx + A\ 2 x 2 A 2 2lXl +Al 2 x 2 

= xjdetiA 1 ) + x 1 x 2 det(A 2 ) + x 2 x x det{A i ) + x\det(A i ) ^ 0. (30) 

Here the new matrices A 1 (i — 1, ... ,4) are constructed as all possible combinations of the columns of the initial 
matrices A k . 

The condition in Eq. (|30[) now reads: if the Jacobian det(^r-) is nonvanishing, then at least one of the determinants 
det Ai (i = 1, ... ,4) is nonvanishing. This statement can be reformulated as a sufficient condition for the degeneracy 
of the transition {x\, x 2 } — » {/i, f 2 }. Namely, if all determinants det A % (i = 1, . . . , A) vanish, than the set of quadratic 
equations H29\) cannot be inverted. 

Now we would like to apply the above criterion to the reaction with two pseudoscalar mesons in the final state. As 
already mentioned, we perform a truncated partial wave analysis, where the amplitude is decomposed over the partial 
wave amplitudes up to some maximum value of the total angular momentum J max . As a set of observables for the 
truncated partial wave analysis it is convenient to choose real and imaginary parts of the coefficients appearing in the 
expansion of the the functions Uj^ pmi, defined in Eq. (I2U|) over the Wigner functions. Rewriting Eq. ([20)1 as 



2 J, 

A'A 

l I'M':IM — 



E E U^ M ^'\^ 2 )Dl, m {R lp ) , (31) 



j— m'ra 

one obtains the observables E/J„v„ /M,A A in terms of bilinear combinations of the partial wave amplitudes t J J£^ 

T I'M'IM,\'\,, . . . i i 7/t/o: , 1\ V* / i\v' I o h I\l2 \ ^ 



C M ' (^.^) = (-i) M ^(2j+i) e w ; J 



J'M'jJMj 



v -v' M'J U -fj,' —M 

J' J j\f JI J 

M'j -Mj m'J I A' - // fi-X 



x tC^,fau 2 )H^fau> 2 ). (32) 
The symmetry property of Eq. (|2T1) leads to the following symmetry of the observables for the interchange A ■<-> A' 

TT I'M'IM,XX' _ ( _i\M' +M+m' +m trjl' ' -M' I-M,X' An* /oo\ 
U jm'm — \ 1 ) \ U j-m'-m ) ■ ^ JJ > 

For the discussion to follow it is convenient to introduce a matrix notation by writing 

n 

uS2 M ' X ' X =Y,^ X ' X) tu, (34) 

k'k=l 

where f3 = (I'M';IM;jm'm) and jfeW = {J^,Mf , z/M,/xW) enumerates the amplitudes. The maximum value n of 
the indices k! and k is equal to the total number of amplitudes t v ^ for a given value of A. The symmetry relation 
in Eq. (|3"51 leads to the matrix relation 



M' ;IM;jm'm,\' X) _ f^\M'+M+m'+m/jj{l'—M';I—M;j—m'—m,XX')\T (35) 

Furthermore, using the property 

E(-ir m (L -i i)=V2l+lS J0 , (36) 
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(/3,A'A) _ x __ x __ x ,_ x , ,.r/r/0»i1\VLl\-lfj[ J J j 

\Mj -Mj 



one obtains for the trace 

fc JMj 

x L > \v -v 0) 2s -n oj \y - (i n - A A - A' 

= (2 Jmax + l^J'O^M'O^/O^MO^jO&n'O^mO^A'A , (37) 

which means that all matrices U^ 1 M ;fM;3 ' m m ' A A ) have a vanishing trace except for the diagonal matrix [/(o°; 00 ; 000 ; AA ). 

The real and imaginary parts of the coefficients U 1 .^^ ' x = /^ ,A A -* + i g^ ,x A ' may now be treated as observables 
for the truncated partial wave analysis. Again we introduce a matrix representation by 



1 n 

= Re U»f™* A ^E ^ » (38) 

1 ™ 

fl Cft*'A) = im ^™ A = J- ^ A ^ , (39) 

fe'fe=i 

where the matrices F^ ,A A ) and G^' A A - ) are respectively the hermitean and antihermitean parts of the matrix U^ ,x A -* 
(symmetric and antisymmetric parts, respectively, in case of a real matrix ) 

F (/3,A'A) = [/(/S.A'A) + 09,A'A) f s ( 40 ) 
G (/3,A'A) _ jj(/9,A'A) _ ;y(/3,A'A) f _ ( 41 ) 

For the application of the criterion of Eq. (|30[) we introduce real and imaginary parts of the amplitudes by 

tj=x j +iy j) j = l,...,n. (42) 

Since an overall phase is arbitrary, we can take one of Xj or yj as zero. For definiteness we set y n = 0, so that 
the amplitude t n is real. Then introducing the 2n — 1 dimension real vector z = {x\, x-2, x n , y\, y2, 2/n-i}i the 
observables /^' A A - ) and g(^' A A ) can be represented by the following real quadratic forms 

2n-l 

a' a) = ^^'^ (43) 

2n-l 

^> A ' A ) = £ , (44) 

»i=i 

where the (2n — 1) X (2n — 1) matrices A^' x A -* and B^' X A ) are determined as 



- (_g (S >* )T o j. («) 

Here, the matrix F^' X A ) is obtained from F^' X A ) by canceling the nth row and the nth column whereas the matrix 
q{P,\ a) j g obtained from G^ ,x A - ) by canceling the nth column. Using the exact expressions of the matrices F/G^' X A ) 
one can easily construct the matrices A^' X A ) and B^' X A ) of Eqs. (|4"5j) and (|4^|) . respectively. 

Now we study in detail only the simplest case, when the partial wave expansion of Eq. (p} is truncated at J max = 1/2. 
Because of angular momentum conservation, requiring |A — fi\ = 1/2 (or /i = A/2), the total number of amplitudes 
^J\u* f° r eacn ^ — is reduced to four. Furthermore, since we exclude from the present consideration linear photon 
polarization, only the coefficients with A = A' appear in the observables in Eqs. (|14l) and (f2"2"j) according to Eqs. ([T 



(fT9"| and (|28|l. Then the subsets of the amplitudes t\l^ 1 corresponding to A = 1 and A = — 1 can be considered 
separately. This is obvious from the fact that in this case the observables are combinations either of the functions 
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TABLE I: Enumeration k of the amplitudes tk = ^J\^ f° r both A-values. The initial nucleon helicity fi is fixed by the angular 
momentum conservation as fi — A/2. 



k 


1 


2 


3 


4 


M 


-1/2 


-1/2 


1/2 


1/2 


V 


-1/2 


1/2 


-1/2 


1/2 



f/gt-P' 11 *) or f / g^'" 1 " 1 ^ ■ Thus we list in Table Q] for both cases the quantum numbers M and v of the amplitudes 
t k = t^/xu^ an( ^ their enumeration k = 1, . . . ,4. 

For the case J m ax = 1/2 one finds ten values for (3 = (I'M'; 00; jm'm) which are listed and enumerated by np from 
one to ten in Table [TT1 One should note that in the absence of target orientation one always has (IM) — (00). The 
corresponding linearly independent 4x4 matrices F^ n ^- xx ^ and G( n P> AA ) are listed in Eqs. (IB1[) and (|B2|) of Appendix IB] 
Those matrices G^ nfi ' ' which are absent in this listing are either zero or depend linearly on the matrices in Eqs. (|B1[) 
and (|B2[) according to the symmetry of Eq. ([35]) . For each A = ±1 the set {F^ ni3 ' XX \ np = 1, .. . , 10} forms a basis in 
the space of symmetric real 4x4 matrices as does the set {G'"' 3:AA - 1 , np = 2, 5, 7, 8, 9, 10} in the space of antisymmetric 
real 4x4 matrices. Except for F^ 1,xx ', all matrices have a vanishing trace. 





TABLE II: Enumeration np of the observables f/g( n P> xx > for /3 = 


(I'M'; 00; jm'm 


and J = 


1/2. 




Tip 


1 2 3 4 5 6 


7 


8 


9 


10 


I'M' 
jm'm 


00 00 00 10 10 10 
000 110 100 000 110 100 


11 

000 


11 
110 


11 
100 


11 
1-10 



Now, in order to find a complete set of observables (up to already mentioned possible discrete ambiguities) , we have 
to construct at least one nonsingular 7x7 matrix using the columns of the matrices A^ 13 '^ of Eq. (|43| and i?( n /3> AA ) 
of Eq. (|46|) . Because of a rather simple form of the constituent matrices F^ nf> ' xy> and G^ nf3 ' XX ' (see Appendix [Bj) . 
it is not difficult to find different combinations of columns which constitute nonsingular matrices. In fact one can 
select almost any set of eight matrices A^f^^ and i?( n /3- AA ). For example, one may take the columns in the following 
combination 

1a(1,A)i4a(2,A)5 3 > i(3 ) a),4a(4,A))3b(2 ) A))4a(5,A)i 1b(5,A)>4a(6,A) , (47) 

where the notation fcA/-B(ns,A) means that one selects the fcth column from the matrix A/ B^ n P' XX > . Now using for the 
differential cross section the expressions in Eqs. (p~5)) and (p~6|) for cr° and E c , respectively, and for the z-component of 
the recoil polarization and in Eqs. (|2"5|) and (124p . respectively, in terms of /W' AA ) and g( n i 3 > xx \ one finds the 
following set of 16 observables (see Table [TTI for the enumeration f / g^ nf3 ' XX ^) 

a°,S c : /( 1 ' AA >, /( 2 ' AA \ /( 3 ' AA \ g^ xx \, A = ±1 , 

pO pc . / (4,AA) j / (5,AA) j / (6,AA) j 5 (5,AA) 5 A = ±1 , (48) 

where the quantities //g("' 9 ' AA ) are determined by the Eqs. ([55]) and (f59"|) . Out of this set one may select any 15 
observables for a complete set. As noted, such a set is only a 'minimal' complete set of observables in the sense, that 
it generally determines the required amplitudes up to possible discrete ambiguities, only. In other words, if one solves 
the corresponding system of 15 bilinear equations, one finds in general more than one solution. 

In order to resolve the remaining ambiguities and thus to find a proper unique solution, additional information on 
other observables is needed. For a proper selection of additional observables we now apply the criterion formulated 
in Ref . jTo| . Given a linear transformation of the amplitudes tk 

tk^t' k =Y,U kk >t k , , (49) 

k'k 

the criterion of Ref. [ici| reads as follows: if there exists a nontrivial transformation U with the property 

rfOU = O (50) 

for all matrices O € {_F ( " 8 ' A A ), G ,< -' lt,:A A '} of the set of selected observables, than for any solution {t k } of this set of 
observables, the amplitudes {t' k } form another solution of the same set, since 

o = X;t;E,ew fc = Y,e k ,(tfou) k , k t k =X)t£<W* , (51) 
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and thus there is a discrete ambiguity. As is mentioned in Ref. [10| this criterion is in general not sufficient since 
it covers only linear transformations — > t' k . Nevertheless, using this criterion one can resolve at least some of the 
possible discrete ambiguities, thus making the general problem easier to solve. 

Since our minimal set includes the observable /( 1 ' AA ) which is proportional to the scalar product 

the transformations U should preserve the moduli of the amplitudes. It is therefore natural to consider primarily 
unitary n x n matrices. The property of Eq. (|50j) is then equivalent to the commutativity of the matrix U with all 
matrices O of the selected set. Application of the criterion in the present case means, that we have to find a nontrivial 
transformation U in the space of unitary 4x4 matrices which commutes with all matrices F^™' 3 ' A A ) and G^ n ^ ,x A - ) of 
the set listed in Eq. (|48|) . 

Such a matrix U is easily found among the diagonal unitary matrices: 




U= , • (53) 



At the same time, it does not commute with any one of the matrices P/G(«£. AA ) for rip = 7, . . . , 10 (see Table HTJ) 
This means that in order to resolve the ambiguity under discussion, the minimal set in Eq. (|48p should be enlarged 
by any of the observables //g( n /3' AA ) belonging to the recoil polarization components P x and P x (or P° and Py). 

It is interesting to note that according to Table U the matrix of Eq. (|53| corresponds to the transformation of 
Eq. (fT2|) which in turn is equivalent to the parity exchange of Eq. (jlip . Therefore, the existence of the ambiguity 
determined by the transformation U in Eq. (|53l) is directly related to our previous conclusion about the necessity of 
measuring the recoil polarization P x or P y in order to separate contributions from states with different parities. In 
order to resolve this ambiguity it is sufficient to enlarge the set of observables in Eq. (l48l) . for example, by /( 8 > AA ). 
However, there exists another ambiguity related to another nontrivial diagonal unitary transformation commuting 
with this enlarged set, namely 



U = 



( e 2 * 3 • • • \ 
1 



(54) 



V • • • ij 



where 023 = 02 — 03 is the relative phase of t-i and £3. 

For the elimination of this last ambiguity one can supplement the existing set by the observable g( 8 ' AA ). It is easy 
to prove that in the case of J m ax = 1/2 the resulting set of observables turns out to be 'fully' complete. To show 
this we firstly note that the magnitudes of all four amplitudes for A = 1 is determined by the set of linear equations, 
corresponding to the four diagonal matrices F^f- 11 ^ with rip — 1,3,4,6 (see Table HI1 and Eq. (|Bip ). Obviously, the 
determination of the absolute squares \ti\ 2 from this set does not involve any discrete ambiguity. Once the magnitudes 
are known, the relative phases 0i3 = 03 — 0i and 024 = 04 — 02 may be unambiguously determined using the four 
matrices F/G {2 ^ and F/G ( - b ' 11 \ This is the only information which may be obtained from the minimal complete 
set. For an unambiguous determination of all four amplitudes we only need one of the remaining phases 023 or 0i4, 
since the second one may always be found from the identity 



013 - 023 + 024 - 014 = . (55) 

As may be seen from Eq. (IB1|) . the relative phase 023 can be extracted from F/G^ 8,11 '. Obviously, the same procedure 

1/2JV. 



can be applied to the subset A = — 1. Thus, in order to unambiguously determine the amplitudes £J\?f J the following 



set of single and double polarization observables is sufficient 

<t°(/o) , E c (/ Q ) , P° Z '%P Z , , P Q ) , P Q X ( C {P X[ , P°) , (56) 

where for Xi one can take either x or y. In Eq. (|56l) we display in parentheses the corresponding notation of Ref. Q 
for the observables. 

One comment with respect to this result is in order. It is clear that all matrix elements i& may be unambiguously 
determined (apart from an overall phase) if the modulus of one amplitude, say, for example, |ii| 2 , and all interference 
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terms t*U are known either directly or through a chain t\tk,t* k ti, . . . .t^U, because such interference terms can be 
expressed as linear combinations of observables. Such a strategy has been discussed and employed in Ref. @. In 
this respect one should note that our set, containing the absolute values of all amplitudes is overdetermined. The 
knowledge of \ti\ 2 for i > 1 is not needed in this case, since these can be obtained from the obvious identity 

W 2 = ^7^, (57) 

for any j. However the structure of our set of equations does not allow the determination of the magnitude of just 
only one partial wave amplitude. 



IV. CONCLUSION 



The present paper is only the first step towards a systematic approach to a model independent partial wave analysis 
of a complete experiment for the photoproduction of two pseudoscalar mesons on a nucleon. The scheme, presented 
here, is based on a model independent formalism of a partial wave expansion developed in Ref. @. 

The procedure for finding a complete set is based on two criteria. The first one, originally developed for deuteron 
photo- and electrodisintegration in Refs. Q allows the elimination of a set of 2n — 1 independent observables. To 
partially resolve possible remaining discrete ambiguities a second criterion from Ref. [Io| is used. In the simplest case 
of truncating the partial wave expansion at J m ax = 1/2 these two criteria turn out to be sufficient for an unambiguous 

determination of the eight amplitudes t^^ 1 '' . The corresponding complete set includes beyond the unpolarized cross 
section, helicity beam asymmetry, as well as recoil nucleon polarization along the z and one of the x or y axes with 
and without circular polarization of the photon beam. It is rather interesting, that the complete set of observables 
for the reactions discussed here necessarily includes recoil polarization in the plane orthogonal to the quantization 
axis. Otherwise, the contributions of the partial waves with the same total angular momentum J but different parity 
cannot be separated. This property distinguishes the present reaction from those with a single pseudoscalar meson 
in the final state, where one can avoid to measure recoil polarization as demonstrated in Refs. [HQ. 

We are aware of the fact that the practical use of the present results for J max = 1/2 is very limited, since waves with 
J = 3/2 appear to be important in both mr and irrj channels even in the low energy region. Furthermore, in the case 
of truncation at J max = 1/2 the matrices F , /G'(™< 3 ' AA ) are very simple and an increase of J mQX to 3/2 will probably 
require not only quantitative but also some qualitative modifications of the approach. Therefore, the generalization 
of this method to higher partial waves will be considered in a forthcoming paper. 

There is, however, an important conclusion, coming from the present study. Namely, whereas a complete experiment 
for a determination of the spin amplitudes of the reactions considered here is quite complicated, because according 
to the analysis of Ref. Q it requires the measurement of a triple polarization observable, the truncated partial wave 
analysis seems to be doable, and thus further developments in this direction may be very promising. 
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Appendix A: Parity separation 

In order to separate the final states of positive parity from those of negative parity corresponding to the parities of 
the intermediate nucleon resonances in the two-step process, we split 0\f J (l p j p m p ) 

O x ^ J {l P 3 p m p ) = 0^ J + (l p j p m p ) + O x A f J -(l p j p m p ) (Al) 
according to the parity (— l) l p +l i of the final partial wave \pq; {il p ^)j p lq)J M) . Explicitly one finds 

0^ J± (l pJp m p ) = t^lj2((-^ lp+lq ±^i L (-^) jp l? P h "L>/2) 



(i i ^ ^{i v \wmp xl \\\) ■ (A2) 



In the same way we split for v = 1/2 the small i-matrix 



,JM _ .JM + , ,JM - / A q\ 

l l/2\fj, — l Xfi ' Xfi > 

where the t J x ™ ± are defined as in eq. © with v = 1/2 and 0\f J± in place of 0\f J . 
For v = —1/2 one obtains 

2 ) = ("D* E (o' 1 4) (m-A A ^)^ mp (-/2)e^-^ 

Ipjp Trip L 

(-l)h+m p0 XLJ {lpjpmph (A4) 

where we have used the symmetry property of the 3j-symbol and the property of the small d-matrices 

O -P) = {-^) 3 - m dl m , m (P) (A5) 



yielding for (3 = tt/2 



From the same property follows 



and thus one obtains 



d\ rn (t/2) = {-lY*- m *dl m (tt/2) . (A6) 
(-l) m °d l ° m M2) = (-lY'd^Ml) , (A7) 



(-D-4,>/2) (* £ _ J M ) - (-l)'^>/2) (* £ / M ) • (A8) 
This leads to the relation 

( _l)^+™ p0 ALJ (/pJpmp) = (_l)A/( ALJ + (Wpmp) _ ALJ- (/pJpmp) ^ _ (Ag) 

Inserting this into eq. (|A4[) leads finally to 

= (-i)^ +M (C + - C") ■ ( A1 °) 

Therefore, the separate parity contributions are given by 

= ± (-1)* +M ^W) ■ (All) 
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Appendix B: Listing of the Matrices F/G (n ^ xx) 



Here we list the ten symmetric and linearly independent matrices {i^(™p< AA ) } rif) — 1, . . . , 10} and the six asymmetric 
matrices {G^™' 3,AA ^ , n@ — 2,5,7,8,9,10} for A = 1. One should consult Table [TTI for the correspondence between 
P = ( I' M'; 00; jm'm) and the enumeration np. 



i?(44i) = I 



1 • • • 

• 1 • • 

■ ■ 1 ■ 

... 1 



f (2,11) = 1 



/-I 



i?(7,ll) 



1 

2y/2 



■ 1 



• 1 
1 • 

1 



ir(8, 



11) 



^(10,11) = _1 



1 • 



The linearly independent matrices G^""f' n ^ are 
/ • • 1 
-1 • • 



G (2,H) 



1 

2s/2 



\ 



■ ■ 1 

-1 . . 



G^ 



id 



G (8,H) = 1 



■ 1 

V- ■ 



• 1 • 

• • 1 



2V2 1 



• 1 • • 

/ . .1 



p(5,ll) = _J_ 
~ Is/2 



-1 



1 

V • -1 



/• • • • 

• • 1 • 

• 1 • • 

v. . . . 



/ . . 1 



1 

2sf2 



■ 1 



V • -1 



• 1 • 



G (9,H) = _L_ 



1 

2V2 



-1 



p (3,11) = 1 



i?(6,H) = I 



-1 • • 

■ 1 • 

■ • 1 



/I 



-1 



■ 1 



r ~ 2V2 



■ -1 



• • -1 



• -1 



£(7,11) _ _J_ 
U 2V2 



-1 



/7(io,n) _ 1 

ur 2 



l ... -I 

V 1 • • 



For A = — 1 and I = M = the matrices are related to the above ones by 

JP I M ' ; 00 ; j m' m, — 1 — 1) ^_^J ^ M f ',00; jm' m, 11) 



(Bl) 



(B2) 



(B3) 



[1] V. Kashevarov, et al, Phys. Rev. C 85, 064610 (2012). 
[2] W. Roberts and T. Oed, Phys. Rev. C 71, 055201 (2005). 

[3] R. L. Workman, Phys. Rev. C 83, 035201 (2011). 

[4] L. Tiator, AIP Conf. Proc. 1432, 162 (2012) |arXiv: 1109.06081 [nucl-th]]. 

[5] V. F. Grushin, in Photoproduction of Pions on Nucleons and Nuclei, edited by A. A. Komar, (Nova Science, New York, 
1989), p. Iff. 

[6] A. S. Omelaenko, Sov. J. Nucl. Phys. 34, 406 (1981). 

[7] A. Fix and H. Arenhovel, Phys. Rev. C 85, 035502 (2012). 

[8] H. Arenhovel, W. Leidemann and E. L. Tomusiak, Nucl. Phys. A 641, 517 (1998). 

[9] H. Arenhovel, W. Leidemann and E. L. Tomusiak, Few Body Syst. 28, 147 (2000) 
[10] W.-T. Chiang, and F. Tabakin, Phys. Rev. C 55, 2054 (1997). 
[11] E. M. Rose, Elementary Theory of Angular Momentum, Wiley New York 1957. 
[12] H. Arenhovel and A. Fix, to be published. 



